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ABSTRACT 

We report on a survey of NH2D towards protostellar cores in low-mass star formation and quiescent 
regions in the Galaxy. Twenty-three out of thirty-two observed sources have significant (S5(t) NH2D 
emission. Ion-molecule chemistry, which preferentially enhances deuterium in molecules above its cos- 
mological value of 1.6 x 10^^ sufficiently explains these abundances. NH2D/NH3 ratios towards Class 
sources yields information about the "fossil remnants" from the era prior to the onset of core collapse 
and star formation. We compare our observations with predictions of gas-phase chemical networks. 

Subject headings: Molecular Processes; ISM: abundances; ISM: general; ISM: clouds; ISM: molecules; 
ISM:deuterium; radio lines: ISM 



1. INTRODUCTION 

Deuterium in molecules is enhanced above the cosmo- 
logical [D]/[H] value of x 10"^ (Linsky et al. 1995) 

by several orders of magnitude in the cold environs of the 
interstellar medium (ISM). Two theoretical pathways can 
account for the abundance of these deuterated molecules. 
Deuterium binds more strongly in molecules than do less 
massive hydrogen atoms, especially at low temperatures, 
in gas-phase reactions between parent species and e.g. Hj^ 
and H2D+ (Dalgarno & Lepp 1984; Watson 1974) and 
other similar molecular ions in dense (n > lO'^cm""^) gas. 
The qualitative understanding that reactions specifically 
with and H2D^ have been enhanced by recent detec- 
tions (Geballe & Oka 1996; Stark, et al. 1999). Their 
abundances are consistent with requirements of gas-phase 
ion-molecule chemistry predictions of most species, in- 
cluding deuterium isotopes. Reactions on grains will also 
contribute to deuterium enhancements. Deuteration has 
a smaller activation energy than does hydrogenation on 
grain surfaces. Additionally the more massive deuterium 
will desorb less quickly from grain surfaces than hydro- 
gen, increasing the interaction time for potential grain- 
induced deuteration (Tielens 1983). Such predictions are 
considered important for the high fractionation 10^^) 
observed for HDCO/H2CO (Loren & Wootten 1985), 
NH2D/NH3 (Walmsley et al. 1987), and CH3OD/CH3OH 
(Mauersberger et al. 1988) towards the hot core Orion- 
KL. In such regions, high temperatures quench the ion- 
molecule deuterium fractionation reactions by rapid de- 
struction of H2D+ and faster reactions (lower energy bar- 
riers) for more abundant species. High observed deuterium 
fractions towards these sources provide a chemical "finger- 
print" of the physical conditions in the pre-protostellar 
gas. Comparisons of deuterated molecules to their more 
abundant isotopes may, then, provide meaningful insight 
into the evolution of the ISM. 



Deuterated ammonia, NH2D, is a useful probe of pro- 
tostellar sources because (1) it can be easily compared to 
NH3 , a well-used tracer of temperature and dense gas con- 
desations in the ISM; (2) it possesses easily observable mil- 
limeter transitions; (3) it has hyperfinc components that 
can be easily resolved and used to estimate column den- 
sity; and (4) its abundance can be compared with chemical 
models which predict variations with temperature, den- 
sity and evolutionary state. NH2D was first observed by 
Turner et al. (1978) in Sgr B2 and Rodriguez-Kuiper et 
al. (1978) in Orion KL. These observations of hot core 
regions yielded ammonia deuterium fractionation levels, 
NH2D/NH3, of order a few times 10~^, significantly higher 
than the local ISM [D]/[H] value of 1.6 x 10"^. Walmsley et 
al. (1987) suggested that longer time-scales for deuterium 
interaction on grain surfaces (versus that of lighter hy- 
drogen atoms) can generate larger deuterated molecular 
abundances. The high temperatures in these cores sub- 
limate or photo-evaporate grain mantles, providing the 
high observed NH2D/NH3 levels. Olberg et al. (1985) 
investigated this phenomenon in the colder, less evolved 
sources L183 and towards the NH3 column density peak in 
TMCl, TMCl NH3, where low temperatures minimize any 
surface contributions and gas-phase chemistry dominates 
molecule formation. They found NH2D/NH3(L183)=0.05, 
but only an upper limit towards TMCl NH3. Tine et 
al. (2000) revisit the difference in ammonia fractiona- 
tion for L183 and TMCl. They find NH2D detections 
towards both sources, with L183 (NII2D/NH3 ^ 0.1) hav- 
ing 5 times the ammonia fractionation level of TMCl 
(NH2D/NII3 = 0.02). Their gas-phase network based 
on the reaction rates of Millar et al. (1991) sufficiently 
explains the abundances and fractionation of ammonia. 
Saito et al. (2000) invoke grain processes to differentiate 
star-filled and star-less ammonia cores with and without 
NII2D detections, respectively. They made this conclu- 
sion despite the fact that Saito et al. (2000) measured the 
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strong fractionation towards L183, with a value similar to 
ours and Tine et al. L183 is a well known starless core 
which lacks any source of radiation to sublimate grains. 
Some of the discrepancy in the interpretation of the ori- 
gin of NH2D may result from the source samples. The 
listed studies have concentrated on either too few sources 
or sources over a small range in physical conditions. 

We present in this paper an experiment to address the 
deuterium enhancements of ammonia, NH2D/NH3, to- 
wards dark cores and Class sources. The former are 
often termed "pre-protosteller" since most do not pos- 
sess strong continuum detectable by instruments such as 
IRAS. In such environs, ion-molecule chemistry fraction- 
ates neutrals and ions such as ammonia and HCO+ (Woot- 
ten 1987) and dominates surface dcuteratated compo- 
nents desorbed back into the ISM. Most of the grain ice 
mantles remain frozen since the ambient temperatures are 
significantly less than 90 K, the sublimation temperature 
of the major ice constituent, H2O. Class sources are 
cold dense cores with spectral energy distributions (SEDs) 
peaked in the sub-millimeter (Td ~ 20 — 50 K) and highly 
coUimated outflows (Andre, Ward-Thompson, & Barsony 
1993). Blake et al. (1995) showed that Class sources such 
as IRAS4 contain a large number of hydrogenated species 
and long-chain molecules that can only be produced by 
grain catalyzation. Thus, Class sources represent an 
ideal location for separating gas-phase effects from grain 
chemistry. The NH2D/NII3 ratio of these sources will re- 
flect what, if any, additional deuteration has occured on 
the grains. 

Our observations of NII2D also reveal that four (and 
possibly six) sources possess asymmetric line profiles. One 
does not expect to see such complex line profiles from rel- 
atively low abundance (and, therefore, low optical depth) 
species such as NH2D. Line-of-sight self-absorption or 
kinematic effects (e.g. multiple cores within the telescope 
beam) can account for such a line profile. In this work, 
we assert that gas-phase chemistry adequately accounts 
for observed deuterium enhancements in ammonia. We 
present a radiative transfer analysis using (multi-layer) mi- 
croturbulent models in §3, and discuss how ion molecule 
reactions can, in fact, effectively explain the observations, 
including self-absorbed profiles in section §5. We discuss 
our main conclusions in section §6. 

2. OBSERVATIONS 



Bester et al. (1983) favor the para-NH2D transition by a 
factor of 3. As in ammonia, one can compute the optical 
depths and column densities using ratios of the hyperfine 
transitions (Mangum, Wootten, & Mundy 1992; Ho & 
Townes 1983). A weaker dipole moment, /i^ = —0.1848 D, 
produces pure rotational transitions within a given inver- 
sion level. It is responsible for the ground state transi- 
tions, J = loi ^ Ooo levels at 332.82251 GHz {v = 0) and 
332.78189 GHz {v = 1). 



2.1. NRAO 12-meter Observations 
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We observed the 85.93GHz Jk^^k 
110.15 GHz Jk^.k, = In 

sitions of NH2D using the NRAO 12 meter^ telescope at 
Kitt Peak, Arizona (Figures 1, 2, and 3). The data are 
from several epochs, mostly dating from 1987 September, 
but a fair portion from 1997 September, using both the fil- 
ter banks and the hybrid digital spectrometer. The recent 
data set was observed with a two-channel SIS junction re- 
ceiver tuned to receive signals in a single side-band. Each 
channel measured orthogonal linear polarizations. We uti- 
lized two banks of filters, 100 kHz and 250 kHz, each 
consisting of 256 channels, resulting in 0.27 km s~^and 
0.68 km s~^per channel resolution. We used the hybrid 
spectrometer in the two-IF mode with 12.5 MHz band- 
width and 512 channels per IF, resulting in a spectral res- 
olution of 24.4 kHz or 0.131km s~^; both linear polariza- 
tions were averaged. All data were obtained in frequency 
switched mode, at a switching interval of 2.5 MHz. The 
excellent quality of the high resolution data from the hy- 
brid spectrometer motivated us to use it exclusively for our 
analysis presented here. Filter bank data were generally 
used to diagnose any systematic problems. 

Sky dip measurements indicated zenith opacities of 0.01 
to 0.4 when few clouds were apparent. The system tem- 
peratures varied between 180 and 300 K for all observa- 
tions reported here. Calibration was achieved by the syn- 
chronous detection of an ambient temperature absorber 
and the sky. The observed line intensities were corrected 
for forward scattering and spillover efficiency rjfss = 0.68 
to place the data on a T^ scale. Column density calcula- 
tions require an additional main beam brightness correc- 
tion of ?7*„ = 0.95. The beamsize of the 12m at 86 GHz 
is 90", and 70" at 110 GHz. Regular pointing checks per- 
formed with planets indicated a positional accuracy better 
than 5". 

2.2. CSO Observations of 

We also report on the detection of the fundamental tran- 
sition of NH2D, J = loi ^ Ooo, at 332.782 GHz (Ortho) 
and 332.822 GHz (Para), towards the dense core L1689N 
with the Caltech Submillimeter Observatory (CSO). The 
data were observed in orthogonal linear polarizations. The 
1024-channel 50MHz acousto-optical spectrometers were 

^ The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative agreement with the National 
Science Foundation. 



2.0.1. The NH2D molecule 

NH2D is a slightly asymmetric top molecule whose ro- 
tational levels are split by inversion doubling, as in its 
more well studied symmetric isotopomcr, NH3, but the 
asymmetry mixes the rotation and inversion. The pri- 
mary dipole moment, /i^ = 1.4652 D, induces rotation- 
inversion transitions which are in turn split by the ^''N 
nuclear quadrupole moment into five hyperfine compo- 
nents (Bester et al. 1983). The inversion -split loi and 
111 levels, corresponding to the /Uc dipole moment, pro- 
duce a symmetric and anti-symmetric transition, one at 
85.926263 GHz (ortho-NH2D; u = ^ 1) and the other 
at 110.153599 GHz (para-NH2D; v = 1 ^ 0) (see Ta- 
ble 1). The relative nuclear statistical weights taken from 
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used. The data were position switched with 30' east or 
west offsets. Tlic bcamwidth at 333 GHz is 20" and the 
main beam efficiency, determined from planet observa- 
tions, is 67% at 345 GHz. 

2.3. NRAO 140-foot Observations 

Observations of the ammonia transitions were made dur- 
ing 1998 June at the 43m telescope of the NRAO in Green 
Bank, West Virginia (Figure 5). The {J,K) = (1,1) 
(23.694506 GHz) and {J,K) = (2,2) (23.722634 GHz) 
transitions were observed to sensitive limits in all of our 
sources. All of the data were observed in orthogonal linear 
polarizations. The 1024 channel Mark IV autocorrelator 
gave a spectral resolution of 0.25 km per channel. The 
data were frequency switched at an interval of 2.5 MHz. 
The variation of gain with elevation was reduced by use 
of the lateral focusing stage. The beamwidth at 23.7 GHz 
is ~ 1' and the beam efficiency is 30%. This beam width 
is similar to that of the 12 meter observations. Thus, dif- 
ferential beam dilution is not considered when comparing 
the NH2D and NH3 column densities, assuming that both 
isotopes emit from similar volumes of gas. Line intensities 
are measured in units of T^*. The forward scattering and 
spillover efficiency, rjiss, is 0.63 (Loren, Evans, and Knapp 
1979). 

3. RESULTS 

3.1. NH3 owdNHzD 

We determined column densities for both NH3 and 
NH2D using a microturbulent radiative transfer model. 
We utilize H2CO as a measure of source physical condi- 
tions. This explictly assumes that H2CO and ammonia 
overlap. In general, H2CO is excited in more dense re- 
gions than NH3. Thus, temperature estimates from both 
species may not be cospatial in a given source. With this 
caveat in mind, we discuss below how physical conditions 
were determined. 

3.1.1. Estimates of Source Physical Conditions 

Chemical models predict significant variations of molec- 
ular abundances with density, temperature, and evolution- 
ary state of the cloud. We model the temperature and 
density structure of each source using NH3 observations 
and formaldehyde observations from Wootten et al. (2001). 
For clarity we discuss some issues involved with H2CO ra- 
diative transfer modeling below. H2CO, an asymmetric 
rotor molecule, is an important probe of molecular cloud 
physical conditions because of: 

• a large number of transitions closely spaced in fre- 
quency but well separated in energy, possessing a 
coupled sensitivity to kinetic temperature and vol- 
ume density when measured through appropriately 
chosen line ratios. 

• a large, uniform, and ubiquitous abundance in the 
ISM (Mangum, Wooten, Wadiak & Loren 1990). 

• multiple transitions through only a few receiving 
systems, minimizing calibration uncertainties. 



One must compare the various transitions to determine if 
an isothermal, iso-density cloud model or one with more 
elaborate physical conditions is necessary to reproduce the 
observed H2CO transitions. 

Several steps arc needed to determine basic physical pa- 
rameters for the sources. The ratio of the NH3 (2,2) and 
(1,1) measurements provides lower limits to the tempera- 
ture of the cold, extended envelopes of each source. We 
use published data with similar resolution to supplement 
some of our (2,2) data with poor signal-to-noise (see col- 
umn 7 of Table 3). Next, we use the formaldehyde data 
which provides sensitivity to the density and the presence 
of any warm gas. A comparison of many H2CO line ra- 
tios using a Large Velocity Gradient model provides initial 
estimates for both density and temperature (Mangum & 
Wootten 1993). Using these initial estimates, we then 
employ a more robust microturbulent radiative transfer 
model (one in which the size scale of velocity turbulence 
is small when compared to the photon mean free path) to 
iteratively calculate the line intensities of the lower opacity 
para-formaldehyde species and obtain the best fit temper- 
ature and density model which reproduces the observed 
brightness temperatures. In some regions, multiple layers 
must contribute to the observed spectrum. 

In general, low excitation transitions of H2CO (e.g. 
Jk-iKi = loi — > loo are compared to estimate the tem- 
perature and density in the cold region. Alternatively, a 
relatively warmer region tends to produce emission from 
higher energy transitions (e.g. Jk-iKi = 5o5 — > 4o4). For 
example, towards IRAS4A, Wootten et al. (2001) find both 
low excitation (/ T^dV[JK_iKi = loi loo] = 4.21 ± 
0.09 K km s~^) and high excitation (/ T^dV[JK_iK^ = 
5o5 ^ 4o4] = 3.63 ± 0.15 K km s"^) lines of H2CO, indi- 
cating that a layer of high density, warm gas coexists with 
colder, low density gas. 

We use these values for n and T as input parameters 
into a microturbulent model for NH2D and NH3. The 
partition function is determined over all energy levels in 
order to eliminate errors due to high temperature approx- 
imations. With the size of the region set by our beam, we 
essentially have one free parameter, the column density, to 
predict the integrated intensities. We match the intensi- 
ties and line profiles of all of the hyperfine components of 
NH2D and NH3 along with estimates of the optical depth 
from the following equation: 

TB(m) ^ 1 - exp[-T{m)] 
Tb(s) 1 - exp[~aT{m)] ^ ' 

to determine the best fit column density. In equation 1 
m and s stand for main and satellite hyperfine lines, re- 
spectively; Tb is the observed brightness temperature; r 
is the optical depth; and a is the ratio of the satellite to 
main hyperfine intensities. We list all results in Tables 
3 and 5 for the deutero-ammonia and ammonia column 
densities. 5(7 upper limits for the line intensities and col- 
umn densities are listed in Table 4. In Tables 6 & 7 we 
list the ammonia fractionation for single-temperature and 
multiple-layer models, respectively. 

The quality of our NH3 data was not sufficient for rea- 
sonable comparison with our NH2D in all cases. In general, 
not all of the hyperfine components were obtained in one 
spectrum because of hardware errors in the velocity offsets. 
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This required us to carefully compare all column density 
estimates based on only a subset of the ammonia hyper- 
fine components, with data of similar spatial and spectral 
resolution found in the literature. We list all references for 
ammonia data in Tabic 5. In all cases wc find good agree- 
ment between our own NH3 column density estimates and 
those available in the literature. This lends confidence to 
our estimates of ammonia deuterium fractionation. 

3.1.2. Para-to-OHho Ratio for NH2D. 

In Table 8 we list the Ortho-to-Para ratio for the 86 and 
110 GHz transitions of NH2D for eight sources. If one com- 
pares the statistical weights, the 110 GHz line strength is 
three times weaker than the 86 GHz transition. Different 
excitation conditions may vary this ratio slightly. Direct 
comparison of the integrated line strengths indicated that 
the ratio J T^dv [86 GHz]/ f T^dviUO GHz] is consistent 
with 3 for all sources we observed in both lines. We also 
determine the column density of ortho-NH2D for compar- 
ison to para-NH2D. Again, we find that the ratios of the 
column densities in the upper levels of both transitions are 
consistent with three. 

3.2. Self absorption 

We find asymmetric spectral features indicative of an in- 
wardly increasing temperature gradient in low mass proto- 
stars in the 86 GHz NH2D transitions towards NGC 1333 
IRAS4A, L1448C, S68 FIRSl, and possibly in NGC 1333 
IRAS4C and S68N. Very high resolution profiles show sim- 
ilar spectral signatures in H2CO, N2H'^, HCO^, and CS 
(Mardones et al. 1997; Gregersen et al. 1997). Line asym- 
metries or skewness were computed using the following 
dimensionless quantity: 

Skewness = ET(V - Vlsr)^AV/(I]TAV) x (2) 

= {(ET(V - Vlsr)'AV)/(STAV)}"'/' 
Here, we take a weighted sum over each of the hyperfine 
components to obtain an average skew of the entire NH2D 
line profile. For a line that possesses most of its emis- 
sion at velocities less than the source's systemic velocity 
("blue- ward" asymmetric), the skewness is negative; it is 
positive for a similar "red-ward" asymmetry. We list our 
calculations for those sources with significantly asymmet- 
ric line profiles in Table 10. This method indicates that 
IRAS4A and L1448C possess the largest blue asymme- 
try The NH2D spectra of IRAS4C and S68N show strong 
blue- ward behavior as well. The relative uncertainty in the 
S68 FIRSl data makes quantifying an asymmetry difficult. 
No sources show significant (Slcr) red- ward asymmetries. 

4. DEUTERIUM CHEMISTRY AND AMMONIA 

Gas-phase synthesis of NH2D is strongly supported 
by observational evidence for Hs"*" (McCall et al. 1998; 

Geballe & Oka 1996) and H2D+ (Stark, et al. 1999) ion- 
molecule chemistry model predictions, and limits to des- 
orption mechanisms in the molecular gas towards our ob- 
served sources. We discuss grain desorption mechanisms 
in §??. For clarity, we briefly discuss the general ion- 
molecule chemistry responsible for NH2D and NH3 pro- 
duction. We use the detailed models of Roberts &: Mil- 
lar (2000) and Turner (2001) to compare our data to 



in §??. We refer the reader to those sources for a more 
extensive discussion. Using the most current set of reac- 
tion rates from the UMIST database (Millar, Farquhar & 
Willacy 1997), Roberts & Millar (2000) determined the 
most important chains for the deuterium injection into the 
ammonia production process. The three general schemes 
considered in the literature are: 

H2 N+ ^ NH+ + H (3) 
H+ N ^ NH+ + H (4) 
H+-hN^NH2+ + H (5) 

Reaction 3 is often rejected for NH3 production because 
experiments performed at very low temperatures indicate 
a reaction endothermicity of 85K (Marquette et al. 1985) 
for thermal conditions. The corresponding reaction with 
HD, on the other had, has an endothermicity of 16K. This 
reaction also has the problem that N"*" is not very abun- 
dant. Herbst, Defrees, & McLean (1987) suggests that a 
nonthermal formation mechanism, where N+ obtains ad- 
ditional translational energy in the reaction 

N2 + He+ ^ (N+)* N He (6) 

can provide the necessary abundance of N+ for equation 
3. A standard set of laboratory-determined reaction rates 
have not been determined. The range of rates lead to a fac- 
tor of 10-100 variation in predicted ammonia abundances. 
Reaction 4 works well only for the formation of NH3, since 
the deuterium equivalent has an energy barrier and limited 
abundance of HD+ . Reaction 5 has received much recent 
debate since its rate was remeasured by Scott et al. (1997) 
to be sufficient for ammonia (and, assuming the same rate, 
NH2D) production. However, this contradicted earlier the- 
oretical work which indicated that double-proton transfer 
was too inefficient (Huntress 1977), even if favored at cold 
temperatures. A more recent measurement for the rate of 
equation 5 indicates that the result of Scott et al. (1997) 
was erroneous. The authors no longer consider it a viable 
channel for ammonia production. 

5. DISCUSSION 

5.1. Comparisons with gas-phase model predictions of 

NH2D. 

We plot in Figure 7 the ammonia fractionation as a func- 
tion of source temperature listed in Table 9. Observed 
values rise to ~ 20 K; sources warmer than this have a 
flat NH2D/NH3 ratio. The range of NH2D/NH3 for mas- 
sive hot cores is included as an arrow on the right side of 
the plot. We include curves (dashed-lines) based on the 
calculations of Roberts & Millar (2000) that bracket a 
range of values for NH2D/NH3. These model points are for 
steady state models at densities of 10^ and 10'''cm^'^. The 
initial conditions of their models include a cosmic D/H and 
HD/H2 ratio of 1.6x10-5 (Linskyet al. 1995), a C/0 ratio 
of 0.4, and a cosmic ray ionization rate of 1.3 x lO^^^s^^. 
The range of model values plotted in Figure 7 includes a 
depletion of C, N, and O abundances by factors of 3 and 6. 
Larger depletions and densities lead to higher NH2D/NH3 
values. 
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The observational result suggests that Class and dark 
core sources have NH2D/NH3 values set by the earlier 
conditions of the parent cloud. Thus, we observe the 
"fossil remnants" of gas-phase chemistry when we study 
NH2D/NH3 towards these sources. The theoretical predic- 
tions of Roberts & Millar (2000), however, appear to con- 
flict with the observed values, even for the coldest sources 
whose formaldehyde models consist of single, isothermal 
zones. The time-dependence of ammonia fractionation 
suggests a way to resolve this discrepancy. Roberts & Mil- 
lar (2000) indicate that NH2D/NH3 values of ~ 0.1 are 
reached very quickly within a few xlO^ yrs, for isother- 
mal 10 K sources at moderate densities of lO^cm^'^ and 
persists for roughly 10^'^ yrs. In a pure gas-phase model, 
this decreases to ~ 10~^ after lO'^"'' yrs, while for a model 
that includes active grain accretion the fractionation rises 
above 10~^. This suggests that all of the sources except 
one (TMCl-CY) arc at least 10"* yrs old. However, such 
age estimates arc unconstrained. Additional deuterium 
fractionation ratios and abundances are necessary in order 
to fully test the validity of such a "chemical chronometer" . 
One example that may help calibrate such a clock is re- 
cent modeling of CO depletion towards L1544 by (CascUi 
et al. 1999). Their dynamical and chemical modeling of 
CO, H^^CO+, and D^^CO+ Hne profiles suggests an age 
of at least lO'' years. They derive a DCO+/HCO+ ratio 
of 0.12 ± 0.02. We find NH2D/NH3 = 0.13 ± 0.02 towards 
L1544 which, according to the model of Roberts & Millar 
(2000), is consistent with the age derived by (Caselli et al. 
1999). 

The two primary Taurus positions, TMCl-CY and 
TMCI-NH3 offers a relative comparison of age estimates 
derived from deuterium fractionation ratios. Pratap et 
al. (1997) mapped the spatial distribution of 34 transi- 
tions of 14 molecules towards the Taurus Molecular Cloud 
ridge, which includes these sources. They argue that 
variations in the emission from carbon-bearing and other 
molecules result from both density and abundance varia- 
tions. Pratap et al. (1997) find especially striking abun- 
dance gradients for SO, HC3N, and CH3CCH, which they 
explain with a small difference in the chemical evolution 
time-scale between the northeast and southwest ends of 
the cloud, or by a small change in the gas-phase C/0 
ratio. Both possibilities effectively reduce the carbon 
abundance, which would explain the lower abundances of 
cyanopolyyenes in the northern regions near the ammo- 
nia peak. Comparisons of NH2D/NH3, DCO+/HCO+, 
DCN/HCN, and N2D+/N2H+ with the chemical model 
of Roberts & Millar (2000) satisfies both of these require- 
ments. However, the age difference is preferred since both 
NH2D and N2H+ are considered to be largely derived from 
N2 . Molecular nitrogen is among the least reactive of neu- 
tral species in the ISM, making it likely to be abundant 
in the dense gas found in the more advanced stages of 
protostellar evolution. 

Two recent studies of NH2D/NH3 by Tine et al. (2000) 
and Saito et al. (2000) support and contradict, respec- 
tively, our conclusion that gas-phase chemistry is re- 
sponsible for the ammonia deuterium ratios of L183 and 
TMCI-NH3. Tine et al. (2000) compare their observa- 
tions with a simple-gas phase network based on the re- 
action rates of Millar et al. (1991). They successfully re- 



produce their observed abundances and deuterium frac- 
tionation ratios. Saito et al. (2000) on the other hand, 
compare NH2D data for several ammonia cores. They in- 
voke grain synthesis of deuterated species even for their 
cold, star-less sources such as L183. Tine et al. (2000) 
observed L134N (1' NW of L183) and TMCl-N (equiva- 
lent to TMCI-NH3). They experienced some ambiguity 
in modeling the spectra with an LVG code and relied in- 
stead on LTE values for determining NH2D column den- 
sities, analysis, and fractionation: N(L134N) = 2 x 10"'^^ 
and N(TMC1 - N) = 1.3 x lO^^; X(L134N) = 1.5 x 10"^ 
and X(TMC1 -N) = 1.2 x 10"^°; (L134N) = 0.18 and 
(TMCl - N) = 0.02. 

We find N(NH2D) equal to 3 x lO^^ ^nd 1 x lO^^ ^^-2^ 
one and zero orders of mag different for LI 83 and TMCl- 
NH3. However, we have beam dilution between effects 
between our 12 m observations and the 30 m survey of 
Tine et al. (2000) and different positions, for L183. Our 
abmidances, columns and fractionations for TMCI-NH3 
are fairly similar. Tine et al. (2000) are able to repro- 
duce their observations with extensive gas-phase modeling 
of the ammonia chemistry, with a few important caveats. 
Moderate to significant (factors of a few to 10) depletion of 
carbon and oxygen (i.e. removal of neutral and atomic de- 
stroyers, and reduction of X^- ) will help achieve the high 
abundances and fractionations (0.02 and 0.1 for L134N and 
TMCl N, respectively). Our results for TMCI-NH3 and 
L183 arc similar, though we find a factor of two smaller 
deuterium fractionation towards LI 83. Additionally, we 
have been able to extend our analysis to other dark cores, 
which also indicates that gas-phase chemistry is most im- 
portant. This suggests that grain processes are unimpor- 
tant for these low mass cores, at least for ammonia forma- 
tion and deuterium fractionation. 

The results of Saito et al. (2000) are discrepant with 
our view that gas-phase synthesis dominates the forma- 
tion of NH2D. They studied NH2D/NH3 for 16 ammonia 
cores from the Benson & Myers (1983) list and found 
that mostly those with IRAS detections show elevated 
or observable NH2D abundances and NH3 fractionation. 
Furthermore, they find larger abimdanccs than expected 
from gas-phase networks for the kinetic temperatures of 
4 sources. They concluded that ammonia is formed and 
deuterated on grains. Several problems exist with their 
analysis. First, a notable exception to their conclusions 
regarding grain chemistry formation of NH2D is L183, 
where they find abundances and fractionations similar to 
our own. Second, they use an LTE estimate of the exci- 
tation of NH2D, which tends to over-predict the column 
density. This is because in an LTE scenario, one requires 
a larger abundance (column density or optical depth of 
NH2D) to achieve the same line intensities than for non- 
LTE methods. Third and last, they did not integrate 
enough on all of their sources, since TMCI-NH3 has a 
clear detection by us and by Tine et al. (2000). Thus, it 
is also possible that for their dark cores without infrared 
sources, the column density of NH2D may be below the 
detection limit of their survey. 

Observational and theoretical comparison of L183 and 
TMCI-NH3 provide important benchmarks for gas-phase 
synthesis of molecules (Swade 1987). L183 and TMCl- 
NH3 possess similar values of n, T, and size, but dis- 
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play very different chemical characteristics. Most chem- 
ical and dynamical models of L183 have concentrated on 
the dearth of cyanopolyyenes and abundance of sulferetted 
molecules. Deuterium fraction may provide additional 
constraints. Ammonia fractionation is larger by a factor 
of two in L183 when compared to TMCI-NH3. Different 
initial abundances between L183 and TMCl have gener- 
ally been ruled out by their relative isolation and similar 
galactocentric distance. TMCI-NH3 should, as a whole, 
be older than L183 based on time evolution models of hy- 
drocarbons; however the same cannot be said for sulphur, 
which is observed to have lower abundances in TMCl re- 
sulting in the conclusion that L183 is older. Swade (1987) 
points out that the conflicting chemical histories of L183 
and TMCl- NH3 can be rectified in terms of larger deple- 
tions in TMCl. In order to achieve this, one requires a 
larger molecular hydrogen volume density, which in cold 
gas leads to more easily depleted material. The slightly 
lower abundance of ammonia (which we observe) is con- 
sistent with a larger depletion in TMCI-NH3. 

5.2. Sources with warm gas and more complex structure 

Our H2CO modeling indicates that several sources 
(L1448IRS3, NGC 1333 IRAS4A, NGC 1333 IRAS2, 
NGC 1333 IRAS7, S68 FIRSl, and S68N) are well fitted 
with two distinct layers with different densities and tem- 
peratures. We do not know for a certainty where all of the 
NH2D and NH3 is located. The reason for this discrepancy 
is the inherent ambiguity in our source models. We are 
able to arbitrarily place the ammonia in cither the warm 
or cold zones, and still suc;cessfully rc;pro(liice the observed 
integrated intensity. Therefore, we know only what the to- 
tal column density and column density ratios are towards 
these sources. The lack of knowledge about the grain con- 
tributions to the observed deuterium fractionation makes 
an accurate census difficult. However, similarity between 
column density ratios among Tables 6 (cold, single layer) 
and 7 (warm, multi-layer) does not suggest an extra warm 
gas contribution is necessary. Those profiles that show 
well-resolved self-reversal provide a tool to break this de- 
generacy. We concentrate on IRAS4A, since it contains 
the most well resolved line profiles, and is thus easiest to 
consider. Higher spectral resolution data are necessary for 
a less ambiguous estimate of the column density distribu- 
tion in the other self-absorbed sources. 

5.2.1. Foreground Absorption in IRAS4A 

An extended, cool envelope may explain the self 
absorbed profile observed towards IRAS4A. We estimate 
a lower limit to the self-absorbing layer of NGC 1333 
IRAS4A. We assume that the NH2D absorbing layer is pri- 
marily cold gas near 25K located along the line-of-sight in 
the central 0.5 km s^^. We estimate the lower limit to the 
absorbing layer column density by assuming that the opti- 
cal depth is 1 and varying the temperature. We find that 
the column density must be at least a few times 10^'' cm^^ 
for temperatures between 10 and 25 K (the temperature 
of the cold layer in our micro-turbulent model). 

This allows us to fix the column density in the outer 
layers of IRAS4 to at least 1 x lO'^"^ cm~^ and estimate 
the total column density with our micro-turbulent model 

^ See the discussion by Mundy & McMuUin (1996). 



by fitting the observed integrated intensity listed in Table 
3. Using this as an initial guess we estimate the column 
density towards IRAS4A in the warm gas is at least a fac- 
tor of 10 less than in the colder gas. The ammonia data 
lacks the spectral resolution to permit a similar analysis. 
If we arbitrarily divide the NH3 column density evenly be- 
tween the warm and cool layers {i.e. assume that NH3 it- 
self is ubiquitous and well-mixed in the ISM), NH2D/NH3 
falls off by about a factor of 10. This is roughly the ex- 
pectation from gas-phase chemical models of deuterium 
fractionation. However, in addition to the lack of spec- 
tral resolution in the NH3 data, we are not resolving any 
spatial structure towards IRAS4A. High spatial resolution 
analysis are necessary to understand the distributions of 
NH2D and NH3. 



5.3. Grain Deuteration 

Grain enhancements may provide an important alter- 
native source for deuterium fractionation of ammonia for 
Class sources. Though we suggest that the flat trend 
in NH2D/NH3 for sources warmer than 20 K indicates 
a gas-phase origin, it is important to discuss the gen- 
eral scheme here. Grain accretion time scales in dense 
cores are of order the collapse time scale, SlQ^yrs. The 
larger mass of deuterium allows it to stick longer to a sur- 
face, and bind more strongly than hydrogen with other 
atoms. Observations of D2CO/H2CO (Turner 1990; 
Castets et al. 1999), HDO/H2O (Jacq ct al. 1990), and 
CH3OD/CH3OH (Mauersberger et al. 1988), underscore 
the importance of grain-based molecular formation and 
fractionation, although more recent gas-phase synthesis 
models (e.g. Roberts & Millar 2000) suggest otherwise. 

A grain origin for NH2D is suggested based on the simi- 
larity of the fractional abundances of NH2D in TMCl ver- 
sus the Orion Hot Core. A review of the abundances using 
standard methods is appropriate at this point. Walmsley 
et al. (1987) find that the NH2D emission is optically thin 
in the Orion-KL region. They obtain N(NH2D)=1.54 ± 
0.4 X 10^^. Using the dust continuum measurements of 
Masson et al. (1985) to estimate the NH2D fractional 
abundance, one finds X[NH2D]=1.01 - 6.21 x lO'^". Us- 
ing the C^^O value from Ungerechts et al. (1997) and the 
conversion factor for carbon monoxide to molecular hy- 
drogen developed by Frerking, Langer & Wilson (1982), 
N(H2) is a factor of 10 smaller than the dust value. There- 
fore, the abundance of deuterated ammonia is equivalently 
larger by ten. Towards TMCl, only C^^O data is avail- 
able, and we find NH2D/H2 = 4.91 ± 1.63 x 10"^°, quite 
similar to the Orion Hot Core. Abundance comparisons, 
in general, are problematic, especially between physically 
different regions such as TMCl and the hot core associated 
with IRc2. CO to molecular hydrogen conversion factors 
and line-of-sight complexity towards a source are difficult 
to properly estimate. Therefore, abundance studies should 
be approached with some caution'^. 

An important distinction between our source list and 
hot cores used in studies such as Walmsley et al. (1987) 
and Jacq et al. (1990) is the temperature. Desorption 
mechanisms will be largely ineffective at removing any de- 
pleted or surface-chemistry products in dark clouds, unless 



7 



low-tempcraturc, non-thermal desorption such as H2 for- 
mation on grains (which releases ~ 4.5 eV of energy) oc- 
curs rapidly. This will only occur for low density regions, 
where H2 formation theoretically prevents the condensa- 
tion of CO and water on small grains and PAHs (Dulcy 
& Williams 1993). Even for the more evolved Class 
objects included here, grain chemistry can be excluded 
because of the sublimation temperature of ices. Grain 
mantles are dominated by water ice (e.g. Ehrenfreund & 
Charnlcy 2000), whose sublimation temperature is 90K. 
Since our sources are all cooler than this, grain sublima- 
tion is not likely. In fact, thermal grain desorption at 
low temperatures and densities above 10"'' cm^"^ proves 
quite incfhcicnt, even in the most optimistic models where 
tunneling on small grains occurs, reaction exothermici- 
ties exceed binding energies, and H/H2 is large enough 
to quickly saturate C, O, and N atoms. We discuss several 
other grain desorption mechanisms below, under the as- 
sumption that surface deuterium enhancements may con- 
tribute to the observed NH2D/NH3. These include: ther- 
mal, cosmic-ray induced grain heating, photo-desorption, 
coUisional (grain-on-grain violence). Impulse heating of 
grains (temporarily increasing the temperature of a sin- 
gle grain above the threshold for thermal desorption) is 
achieved with X-rays and cosmic rays. Leger et al. (1985) 
investigated this affect in regions of varying density and 
visual extinction. For n < 10^ cm~^ and Ay < 5, desorp- 
tion cfhcicntly removes CO and even water ammonia ice 
mantles quite readily. However, for more dense regions, 
the water-ammonia lattice will remain, since for AySfew 
UV photons cannot penetrate them. Wootten, Loren, & 
Snell (1982); Butner, Lada, & Loren (1995); Wilhams 
et al. (1998) have shown that X-ray contributions appear 
to be minimal since the electron fraction measured via 
DCO+/HCO+ varies little for sources with and without 
stars. Therefore, our sources would have little cosmic ray 
desorbed input from abundances of ammonia either stored 
or produced on grains. 

6. CONCLUSIONS 

We observed NII2D hyperfine transitions towards 32 
protostellar and pre-stellar sources with similar spatial 
resolution. We find nearly a 70% detection rate for the 
NII2D lines at 85.9 GHz. The observed abundances of 



NH2D when compared to single dish NH3 observations 
of similar beam-width, indicate that the deuterium frac- 
tionation is large, 10-3<]s^jj2d/NH3;£10-i. The observed 
NH2D/NH3 values generally exceed or equal those seen in 
hot core regions as well as in warm, embedded condensa- 
tions in otherwise low luminosity sources with large, cold 
envelopes. Sources with Tr-~20 K follow gas-phase pre- 
dictions for NH2D/NH3; for TfiT > 20 K the trend flattens, 
indicating that Class sources have not begun to destroy 
NH2D. NH2D/NH3 ratios reflect the "fossil" remnants of 
gas-phase synthesis. We additionally conclude that grain 
formation of deuterated ammonia is not necessary to ex- 
plain NH2D/NH3. 

Some evidence exists that dynamical ages of protostel- 
lar sources can be derived from comparisons of deuterium 
fractionation and chemical models of star forming cores. 
However, this notion is currently rather speculative. We 
are able to address relative evolutionary differences be- 
tween sources such as TMCl-CY, TMCI-NH3, and L183 
successfully. 

Self-reversed profiles indicate that the overall story for 
NH2D is far more complicated, however. Indeed, future 
careful modeling and observations will attempt to address 
the coupled influences of chemistry and dynamical evolu- 
tion in IRAS4. Nonetheless, NH2D provides important 
insights into deuteration in the ISM. Deuterium fraction- 
ation not only varies among sources of different physical 
attributes, but also among diff^crent molecules. Thus, a 
full census of deuterium fractionation seems necessary in 
order to test astrochemical models. Furthermore, with 
deuterated molecules now providing information on the 
collapse regions of protostars, new and potentially pow- 
erful tools for examining the history of core collapse via 
chemical evolution models can provide checks on more tra- 
ditional methods. 
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versity of Illinois. 
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Fig. 1.— 85.9 GHz NH2D spectra. The spectral resolution is 0.13 km s-^/channel except for NGC1333 DCO+ and TMCl NH3 where it is 
0.35 km/s. A linear baseline is removed from the data, revealing the hyperfine components. 
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Fig. 2. — Same as figure 1. Here TMCl CY, L183, L1689N and L1551 are observed with the coarser 0.35 km/s resolution. 
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Fig. 4. — 110.1 GHz NH2D spectra. The spectral resolution is 0.35 km s ^/channel. A linear baseline is removed from the data, revealing 
the hyperfine components. 
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Fig. 6. — Plot of the 333 GHz J = In — » Oqo fundamental transition NH2D observed with the Caltech Submillimeter Observatory. The 
spectral resolution is 0.04 km for this 1.5 hour spectrum. Both sets of lines near 332.820 and 332.782 GHz are heavily blended. The 
integrated intensity is consistent with an LTE distribution. 
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Fig. 7. — Plot of ammonia fractionation versus temperature for sources with single temperature H2CO models (crosses). We include the 
predicted ammonia fractionation from the steady-state gas phase models of Roberts & Miller (2000) as upper and lower ranges marked with 
dashed lines. Additionally, we plot the range of values found for hot cores by various authors {N.B. the temperature choice is arbitrary) 
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0.417 


110.153599 


J 


= 111 


^ loi 


1-1 


0.083 


110.153599 








1-2 


0.139 


110.154222 








1-0 


0.111 


110.155053 








0-1 


0.111 


332.8215595 


J 


= loi 


^Ooo 


2-1 


0.556 


332.8224149 








1-1 


0.333 


332.8229853 



tTine et al. 2000; Townes & Schawlow 1955 
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Table 2 
Source List 



Source il. A. (1950) DEC (1950) Vlsr 

h m s ° ' " kms-i 



L1448 IRS2 


03 


22 


17.9 


30 34 41 


4.0 


L1448 NW 


03 


22 


31.1 


30 35 14 


5.4 


L1448 IRS3 


03 


22 


31.9 


30 34 45 


4.7 


L1448 C 


03 


22 


34.4 


30 33 35 


5.4 


NGC 1333 DCO+ 


03 


26 


3.60 


31 04 42 


7.0 


NGC 1333 IRAS7 


03 


26 


6.90 


31 08 28 


7.0 


NGC 1333 IRAS2 


03 


25 


52.6 


31 04 30 


7.0 


NGC 1333 IRAS4A 


03 


26 


4.78 


31 03 14 


7.0 


NGC 1333 IRAS4B 


03 


26 


7.00 


31 02 52 


7.0 


NGC 1333 IRAS4C 


03 


26 


8.10 


31 03 37 


7.0 


IRAS 03282 


03 


28 


15.2 


30 35 14 


7.0 


IRAS 03367-IC 348 


03 


36 


47.1 


31 47 29 


7.0 


L1551 


04 


28 


40.2 


18 01 42 


6.5 


TMCl AM 


04 


38 


19.5 


25 42 29 


6.0 


TMCl CY 


04 


38 


38.0 


25 35 45 


6.0 




UO 


UU 


04.0 


Qo /in nn 
oZ 4U UU 


v o 
1 .y 


L1544 


05 


01 


15.0 


25 07 00 


7.1 


OMC2 


05 


32 


58.0 


-05 12 11 


11.0 


HHl 


05 


33 


52.0 


-06 47 09 


10.0 


NGC 2264 


06 


38 


24.9 


09 32 29 


8.0 


L183 


15 


51 


30.0 


-02 43 31 


3.0 


VLA1623 


16 


23 


25.0 


-24 17 47 


4.0 


OPHBl 


16 


24 


09.0 


-24 22 49 


3.0 


OPHB2 


16 


24 


26.3 


-24 19 49 


3.0 


IRAS 16293 


16 


29 


21.0 


-24 22 16 


4.1 


L1689N 


16 


29 


27.6 


-24 22 08 


4.0 


L63 


16 


47 


17.0 


-18 00 00 


6.0 


S68FIR 


18 


27 


17.5 


01 13 23 


8.0 


S68N 


18 


27 


15.9 


01 14 49 


9.0 



Table 3 

Deutero-Ammonia Microturbulent Column Density Fits 



Source 


AV 






n(H2)t 


N(NH2D) 


Temp. 




km 


K km s^^ 


K 


10"^ cm"3 


10^^ cm" 2 


Ref. 


L1448 IRS2 


0.49 


0.84 ±0.07 


20 


0.1 


7.60 ±0.61 


OTingcr et al. 1999 


L1448 NW 


1.38 


3.01 ±0.60 


30 


0.1 


36.0 ±7.20 


Barsony et al. 1998 


L1448 IRS3 


1.10 


1.39 ±0.09 


20-50 


0.1-5 


8.14 ±0.55 


This paper 


L1448 C 


0.93 


0.81 ±0.26 


40 ±10 


1 


9.30 ±2.96 


This paper 


NGC 1333 DCO+ 


2.60 


3.77 ±0.96 


20 


0.6 


45.0 ± 11.5 


This paper 


NGC 1333 IRAS2 


1.31 


0.46 ± 0.09 


20-85 


0.4-1.26 


2.62 ±0.49 


Wootten et al. 1999 


NGC 1333 IRAS7 


2.31 


2.30 ± 0.46 


40 


0.3 


2.60 ±0.52 


Lcfloch et al. 1998 


NGC 1333 IRAS4A111 


1.50 


3.98 ± 1.19 


25-50 


0.7-2 


21.8 ±6.50 


Wootten et al. 1999 


NGC 1333 IRAS4C 


1.10 


1.87 ± 0.11 


15 


0.1 


15.4 ±0.90 


This paper 


IRAS 03282 


0.90 


0.76 ± 0.23 


20 


0.1 


7.00 ±2.12 


Wootten et al. 1999 


L1551 


1.35 


1.04 ±0.11 


50 


0.1 


22.5 ±2.40 


Wootten et al. 1999 


TMCl AM 


1.12 


0.44 ±0.12 


10 


0.1 


9.77 ± 2.57 


This paper 


L1544 


0.96 


2.76 ±0.12 


10 


0.3 


26.0 ± 1.10 


This paper 


OMC2 


1.66 


2.55 ± 0.37 


20 


0.1 


220.0 ±32.0 


Batrla et al. 1983 


HHl 


1.68 


1.04 ±0.33 


15 


0.1 


100.0 ±32.0 


Martin-Pint. (1987) 


NGC 2264 


3.10 


0.41 ±0.08 


25 


0.1 


34.4 ± 5.20 


dc Bois. et al. 1996 


L183 


0.30P1 


2.40 ± 0.42 


10 


0.01 


30.0 ± 11.1 


This paper 


VLA 1623 


1.04 


1.92 ±0.17 


20 


0.1 


17.0 ± 1.50 


AWB93 


Oph B2 


1.33 


1.90 ±0.40 


13 


0.1 


9.00 ± 1.89 


Lefloch et al. 1998 


IRAS 16293 


0.81 


1.24 ±0.12 


30 


0.6 


4.50 ± 0.30 


van Dish, ct al. 1995 * 


L1689N 


1.22 


6.46 ± 0.45 


8 


0.1 


160.0 ± 11.0 


Wootten (pc) 


S68 FIR 


1.58 


1.23 ±0.43 


20-100 


0.3-1.26 


34.0 ± 12.00 


McMuUin et al. 1999 


S68N 


1.46 


1.32 ±0.15 


15-75 


0.3-2 


7.82 ± 0.87 


McMullin et al. 1999 



t Sources with two values list the H2CO fits for density and temperature in two layer models 
tProm the fit to their cold, extended layer 

[^iThe error listed contains an estimate to the error in the linewidth of 0.3 km s~^. 
PI Using line fits from high resolution 30 kHz spectrum. 
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Table 4 
NH2D 5t7 Upper Limits 



Source 




N(NH2D) 




K km s ^ 


cm " 


IC 348 


0.094 


2 X 10^^ 


L1512 


0.024 


1 X 10^2 


TMCl CY 


0.007 


5 X 10" 


Oph Bl 


0.045 


1 X 10^2 


L63 


0.060 


1.5 X 10^2 



iN.B. Upper limits from Tk = 20 K 
and n(H2) =1 x 10^ 



Table 5 

Ammonia Microturbulent Column Density Fits 



Sourcet 


AV 
km 


K km s"^ 


K 


n(H2)t 
lO'' 


N(NH3) 
1014 cm-2 


Ref. 


L1448 IRS2 


0.87 


6.21 ± 1.32 


20 


1 


3.23 ± 0.69 


This paper 


L1448 NW 










2.00 ± 0.20 


BaJTsony et al. 1998 


L1448 IRS3 


1.22 


14.74 ± 0.74 


20-50 (30) 


0.1-5 


4.53 ± 0.23 


This paper 


L1448 C 


1.37 


11.07 ±2.35 


40 ± 10 


1 


3.23 ± 0.69 


This paper 


NGC 1333 DCO+ 


1.56 


15.07 ± 1.56 


< 30 (38) 


0.6 


3.71 ± 0.38 


This paper 


NGC 1333 IRAS4A 


1.52 


12.09 ±1.61 


25-50 


0.7-2 


3.08 ± 0.41 


Blake et al. 1995 


NGC 1333 IRAS4C 


1.44 


9.73 ± 3.50 


25 (40) 


0.1 


2.30 ± 0.82 


This paper 


IRAS 03282t 








10.0 ±2.0 


BMP91 


IC 348t 










2.50 ± 1.50 


BGK87 


L1551 


1.27 


0.78 


5.33 


0.1 


2.60 ± 0.26 


This paper 


TMCl NH3 


0.81 


5.48 ± 1.37 


20 


0.1 


1.55 ± 0.39 


Pratap et al. 1997 


TMCl CY 


0.72 


3.34 ± 0.45 


20 


0.1 


0.91 ±0.12 


Pratap et al. 1997 


L1512t 










7.94 


BM83 


L1544 


0.73 


5.04 ± 0.89 


10 


0.3 


1.97 ±0.35 


This paper 


OMC2t 










75.9 ±8.0 


CW94 


HHlt 










5.00 ± 0.38 


MP87 


NGC 2264t 










8.00 


Kriigcl ct al. 1996 


L183 


0.73 


7.52 ± 1.45 


10 


0.01 


5.19 ±0.91 


Olbcrg ct al. 1983 


VLA1623t 










2.00 


Woottcn et al. 1999 


Oph Bit 






19 




3.45 ± 1.75 


MP83,ZBW84 


Oph B2 


1.20 


4.08 


13 


0.1 


0.96 


MP83 


IRAS 16293t 










15.0 ±5.0 


Mundy ct al. 1995 


L1689N 


1.00 


0.55 


8 


1 


18.0 


Woottcn (pc) 


L63t 










10.0 


BM83 


S68 FIR 


1.25 


8.83 ± 2.33 


20-100 


0.3-1.3 


2.71 ±0.72 


This paper 


S68N 


1.24 


5.17 ±0.77 


15-75 


0.3-2 


1.45 ± 0.22 


This paper 



tNHa column densities for Sources without listed line parameters axe taJjen from the noted references, 
t Sources with two values list the density and temperature in two layer models 

lMP83=Martin-Pintado, Wilson, Gardner, & Henkel (1983) 
2ZBW84=Zcng, Batrla, & Wilson (1984) 
^MP87=Martin-Pintado & Cernicharo (1987) 
4BM83=Benson & Myers (1983) 
5CW94=Cesaroni & Wilson (1994) 
'3BGK87=Bachiller, Guillateau, & Kahane (1987) 
'^BMP91=Bachiller, Martin-Pintado, & Planesas (1991) 
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Table 6 

Ammonia Fractionation for Single Temperature and density models 



Source 


Tk 


n(H2) 


NH2D/NH3 




K 







L1448 IRS2 


20 


1 


X 


10*5 


0.024 ±0.005 


L1448 NW 


30 


1 


X 


10*' 


0.09 ±0.02 


L1448 C 


4Q+30 


1 


X 


106 


0.029 ±0.006 


NGC 1333 DCO+ 


< 30 


1 


X 


iqb 


0.10 ±0.03 


NGC 1333 IRAS4C 


25 


1 


X 


10^ 


0.067 ±0.013 


IRAS 03282 


20 


1 


X 


iqs 


0.007 ±0.003 


L1551 


50 


7 


X 


iqs 


0.087 ±0.013 


TMCl NH3 


20 


1 


X 


iqs 


0.06 ±0.02 


LI 544 


10 


1 


X 


10''' 


0.13 ±0.02 


OMC 22 


20 


1 


X 


10'' 


0.029 ±0.006 


HHl 


15 


1 


X 


10-' 


0.074 ±0.024 


NGC 2264^ 


25 


1 


X 


10^5 


0.043 ± 0.007 


L183 


10 


1 


X 


10* 


0.058 ± 0.010 


VLA 1623 


20 


1 


X 


106 


0.085 ± 0.008 


Oph B23 


13 


1 


X 


iqb 


0.090 ± 0.019 


IRAS 16293 


30 


6 


X 


iqb 


0.0028 ± 0.001 


L1689N 


8 


1 


X 


iqb 


0.089 ± 0.006 



INH3 data from Kriigcl ct al. (1996). 
^NHa data from Cesaroni & Wilson (1994) 
^Martin-Pintado et al. 1983 



Table 7 

Multi-layer Models of Ammonia Deuteration 



Source 


n 


T 


NH2D/NH3(i) 




cm~^ 


K 




NGC 1333 IRAS4A 


7 X 10-"^ 


25 


0.07 




2 X 10*' 


50 




S68 FIR 


3 X lO'' 


20 


0.13 




1.26 X 10*' 


100 




S68N 


3 X lO'' 


15 


0.05 




2 X 10*' 


75 




L1448 IRS3 


1 X 10^ 


20 


0.02 




5 X 10« 


50 





Summed over both layers. 



Table 8 
NH2D Ortho-Para Ratio 



Source 


/ TAdv{llQ GHz) 
K km 


/ TAdv{86 GHz) 
K km 


/ TAdv(86GHz) 
J TAdvillOGHz) 


N(86GHz) 
N(llOGHz) 


N1333 


0.598 ±0.142 


2.447 ± 0.623 


4.09 ± 1.43 


4.11 ± 1.30 


LI 544 


0.295 ±0.144 


1.792 ± 0.075 


6.08 ± 2.98 


6.10 ±3.02 


OMC2 


0.90 ±0.13 


1.66 ± 0.24 


1.84 ±0.38 


1.84 ±0.90 


HHl 


0.41 ±0.15 


0.68 ± 0.21 


1.66 ± 0.79 


1.66 ± 0.92 


L183 


0.40 ± 0.06 


1.56 ± 0.27 


3.95 ± 0.92 


3.96 ± 0.92 


OPHB2 


0.42 ±0.11 


1.23 ± 0.26 


2.94 ± 1.01 


2.95 ± 1.20 


L1689N 


1.. 59 ±0.21 


4.19 ± 0.29 


2.63 ± 0.39 


2.61 ± 0.45 


LI. 1.^,1 


(),:-!() ± 0,21 


0,()l) ± 0,07 


2,31 ± 1,85 


2,31 ± 1,80 



Table 9 

Ammonia Deuteratium Fractionation 



Source 


Tk[K] 


NH2D/NH3 


DCO+/HCO+t 


DCN/HCNt 


L1448 IRS2 


20 


0.024 ± 0.005 






L1448 NW 


30 


0.09 ± 0.02 






L1448 IRS3 


20 


0.018 ± 0.001 






L1448 C 


40 


0.029 ± 0.011 






NGC 1333 DCO+ 


30 


0.10 ±0.03 


O.O24I 




NGC 1333 IRAS4A 


25-50 


0.071 ± 0.023 






NGC 1333 IRAS4C 


15 


0.067 ± 0.024 






IRAS 03282 


20 


0.007 ± 0.003 






IC 348 


20 


< 0.008 






L1551 


50 


0.087 ±0.013 


O.O35I 


0.016 ±0.001^ 


TMCl AM 


10 


0.06 ± 0.02 


0.0272 


0.0223 


TMCl CY 


10 


< 0.006 


0.00004 




L1512 


20 


< 0.006 






L1544 


10 


0.13 ± 0.02 


0.12 ±0.02 




0MC2 


20 


0.029 ± 0.006 




O.Ol^ 


HHl 


15 


0.074 ± 0.024 






NGC 2264 


25 


0.043 ± 0.007 


O.OI7I 




L183 


10 


0.058 ± 0.010 


O.O7I 




VLA1623 


20 


0.085 ± 0.008 






Oph Bl 


19 


< 0.003 






OPH B2 


13 


0.09 ± 0.02 




0.013 


IRAS 16293 


30 


0.003 ± 0.001 






L1689N 


8 


0.089 ± 0.006 






L63 


20 


< 0.002 






S68 FIR 


20-100 


0.126 ±0.055 


0.005 


0.010 


S68N 


15-75 


0.054 ± 0.010 




0.009 



tAssumes i^c/isc = 60 
^Williams et al. 1998 
^Butner et al. 1995 
^Greason 1986 



Table 10 
Observed Asymmetries 



Sourco 


NH2D''' 


Skcwncss 


CS 


H2CO 


HCO+ 


HCO+ 




■/ = In - loi 




,/ = 2 - 


■/ = 2i2 - lii'' 


./ = :i — 2' 


./ = 1 — 


L1448C 


Blue 


-0.014 ±0.004 


Blue 


Red 


Red 


Red 


IRAS4A 


Blue 


-0.38 ±0.09 


Blue 


Blue 


Blue 


Blue 


IRAS4C 


Blue 


-0.010 ± 0.004 










L1551 IRS5 


None 




None 


Red 


None 


None 


S68N 


Blue 


-0.704 ± 0.194 


None 


Blue 


None 


None 


S68 FIR 


Red?"* 




Blue 


Red 


Red 


Blue 



''This paper; only the 86 GHz data displays these asymmetries 
•"Mardones et al. (1997) 
^Gregersen et al. (1997) 

''This is a fit by eye, not by a strict Gaussian fit to the line features 



